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Clinical Article 



Clinical Application of Gamma Knife Dose Verification 
Method in Multiple Brain Tumors : Modified Variable 
Ellipsoid Modeling Technique 
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Objective : The Leksell Gamma Knife® (LGK) is based on a single-fraction high dose treatment strategy. Therefore, independent verification of the 
Leksell GammaPlan® (LGP) is important for ensuring patient safety and minimizing the risk of treatment errors. Although several verification tech- 
niques have been previously developed and reported, no method has ever been tested statistically on multiple LGK target treatments. The purpose 
of this study was to perform and to evaluate the accuracy of a verification method (modified variable ellipsoid modeling technique, MVEMT) for mul- 
tiple target treatments. 

Methods : A total of 500 locations in 1 0 consecutive patients with multiple brain tumor targets were included in this study. We compared the data 
from an LGP planning system and MVEMT in terms of dose at random points, maximal dose points, and target volumes. All data was analyzed by t- 
test and the Bland-Altman plot, which are statistical methods used to compare two different measurement techniques. 

Results : No statistical difference in dose at the 500 random points was observed between LGP and MVEMT. Differences in maximal dose ranged 
from -2.4% to 6.1 %. An average distance of 1 .6 mm between the maximal dose points was observed when comparing the two methods. 
Conclusion : Statistical analyses demonstrated that MVEMT was in excellent agreement with LGP when planning for radiosurgery involving multiple 
target treatments, MVEMT is a useful, independent tool for planning multiple target treatment that provides statistically identical data to that pro- 
duced by LGP. Findings from the present study indicate that MVEMT can be used as a reference dose verification system for multiple tumors. 
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INTRODUCTION 

Leksell Gamma Knife® (LGK) system has been widely used 
to treat various brain diseases as a low-morbidity therapeutic 
modality 1 ' 2 * 6 ' 12 ' 13) . Since, this is based on a single-fraction high 
dose treatment strategy, independent verification of the Leksell 
GammaPlan® (LGP) is important for ensuring patient safety 
and minimizing the risk of treatment errors 10,13 '. Previous meth- 
ods for verifying LGP by various algorithms were reported by 
Tsai et al. 13) , Marcu et al. n) , Hur et al. 9) and Zhang et al. 14) . 

However, these verification methods were available only un- 
der a limited set of conditions. The method by Zhang et al. 14) 
could calculate only treatment times, which were accurate only 
without plugs. Other methodologies only tested around the 



center of the target 13 '. Additionally, these studies uniformly 
lacked statistical analysis. Instead of showing statistical agree- 
ment between their method and LGP findings, only the size of 
errors was reported 11 ' 13,14 '. Another significant limitation was 
that previous studies were tested exclusively on single treatment 
target. None of the previous studies showed accuracy for multi- 
ple targets. This means that findings from these previous inves- 
tigations are not applicable to all LGK procedures, because a 
large portion of LGK treatment involves multiple targets. 

In the present study, we have attempted to overcome these 
limitations. To do this, we determined for the first time whether 
a verification method could be applied in multiple target treat- 
ment accurately. Multiple target treatment refers to treatment 
plan for multiple brain lesions, such as multiple brain tumors. 
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Additionally, the accuracy of the verification method was val- 
idated by an expert in statistical analysis, which also was not 
tried so far. 

MATERIALS AND METHODS 

The modified variable ellipsoid modeling technique (MVEMT) 
was used as a verification method for our study 9 '. For this, the pa- 
tient's skull is viewed as an ellipsoid with the center at the mam- 
millary body. The ellipsoid is expressed as following equation 

[(x s -l) cos (3+(y s -m) sin (5 cos cc+(z s -n) sin |3 sin a] 2 

[-(x s -l) sin (3+(y s -m) cos |3 cos a+(z s -n) cos |3 sin a] 2 

+ P- 

[-(y s -m) sin a+(z s -n) cos a] 2 
+ j 1 

c 2 

In which, (x s , y s , z s ) are the coordinates of the shot center co- 
ordinate system (SCCS), (1, m, n) are the coordinates of the 
mammillary body in the SCCS and a=y-90°, |3 is the angle ro- 
tated counterclockwise about the z' s -axis of the skull. A single 
beamlet from the 201 collimators was defined as follows 



y P j hi ,. . 

y = — x,z= — x (1=1, 
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,201) 



where (x pJ , y Pii , Zpj), (i=l, •-, 201) are the geometrical locations 
of the 201 collimators in the SCCS. Dose rate for a single beam- 
let at the arbitrary point (p) is determined by the equation 

D call8 (0)-uve ln2t/, '° / d sf \ 2 A 
D^p; Wc ; t)= ' (^J -e HW-wo). PW , F(a)c . r . } 

in which, D caUS (0) is the calibration dose rate for the 18 mm 
collimator helmet, co c is the collimator factor, d; is the "i" beam- 
let depth, PWi is the plug weight (either 0 or 1), and F(co c ; rO is the 
dose profile function related to dose distribution on the trans- 
verse plane of the beam axis. Therefore, the dose rate for one 
shot at point (p) is calculated with the formula 

201 

D(p;co c ;t)=^D 1 (p;a) c ;t) 

i=l 

and, the dose at (p) from multiple shots (n) is expressed by 
the formula where tj (LGP) is the shot time calculated by LGP. 

n 

D(p)=£D s (p; u) c ;tK(LGP) 

s=l 

The equations were programmed using the hypertext prepro- 
cessor PHP scripting language (http://www.php.net). Calcula- 
tion for all matrix points and the corresponding data were 
logged using MySQL® software (http://www.mysql.com) and 
later queried for necessary parameters. Starting in June 2010, 
we consecutively selected 10 patients with metastatic brain tu- 
mors with 2 to 4 targets. We treated all patients using an LGK 
stereotactic radiosurgery system. 

Dose for 50 points in each patient 

On the computer screen, we randomly selected 50 points for 
each patient, 40 within the target and 10 were out of the target. 
Those 50 coordinates of matrix along with the dose calculation 



at each point were assessed, by the LGP system. We then en- 
tered the same series of coordinates into the MVEMT program 
that calculated the dose at each location along with maximal 
dose (location) and target volume. All data were subsequently 
compared to the LGP results. 

Maximal dose points 

The maximal dose point refers to the point with the highest 
radiation dose for each round of Gamma Knife® treatment. We 
calculated the dose and dose rate at each of the 31x31x31 or 
29791 point for one matrix. Dose rates were calculated while 
accounting for the contributions from the 201 beamlets. Since 
all multiple targets were irradiated with each treatment, the 
dose contribution from each shot was also included. 

Target volume 

When a new treatment is being planned, a three-dimension- 
al conformal radiation therapy should be carefully designed 
which accurately conforms the isosurface of a given radiation 
dose to the anatomical boundaries of the tumor in its entire 
three-dimensional configuration. When the maximal dose is 
determined, the number of matrix points whose dose exceed- 
ed 50% of the maximal dose is queried in order to obtain the 
50% dose volume. The target volume is defined as 50% dose 
volume. 

Statistical analysis 

To analyze the agreement of data produced by the two sys- 
tems, LGP and MVEMT, a paired t-test was used (SPSS, ver- 
sion 12.0 : SPSS Inc., Chicago, IL, USA) : ap value of less than 
0.05 was considered statistically significant. We also created a 
Bland- Altaian plot to graphically compare the dose rates for 50 
points as calculated by the two systems 8 '. The Bland-Altaian 
plot is constructed with the MedCalc software program (Med- 
Calc®, Mariakerke, Belgium, http://www.medcalc.be/). 

RESULTS 

All pooled data included that from 500 points (50 points from 
each of the 10 patients). The average dose calculated by LGP 
showed 26.019±8.560 Gy, while the average dose calculated by 
MVEMT revealed 26.075±8.808 Gy. 

Bland- Altaian plots for all individual cases are shown in Fig. 1. 
For all 500 points, the mean difference between doses calculated 
by two methods, LGP and MVEMT, was -0.057 (95% CI) and 
standard deviation was 1.3353. The lower limit (2SD) of differ- 
ence was -2.674 and the upper limit (2SD) of difference was 
2.561. The case 6 revealed a mean difference of 0.692±2.542, de- 
viation of data was much greater than the other cases, and had 
higher magnitude limit values (-4.290 for the lower limit and 
5.674 for the upper limit). When we removed data for the case 
6 from our analysis, better agreement between the two methods 
was observed in which the mean was -0.140 with a standard 
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Fig. 1. Bland-Altman Plots for each individual case. The middle line represents the mean difference between the LGP and MVEMT. The two extreme 
lines are the +1 .96 and -1 .96 standard deviations of the difference between the two methods. MVEMT : modified variable ellipsoid modeling tech- 
nique, LGP : Leksell GammaPlan. 
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Fig. 2. A : Bland-Altman Plot for all pooled data (n=500). The mean difference between doses calculated by two methods, LGP and MVEMT (95% CI) 
was -0.057. B : Bland-Altman Plot for data excluding case 6 (n=450). The mean difference between doses calculated by two methods (95% CI) was 
-0.140. The middle line is the mean difference between LGP and MVEMT. The two extreme lines are the +1 .96 and -1 .96 standard deviations of the 
difference between two methods. MVEMT : modified variable ellipsoid modeling technique, LGP : Leksell GammaPlan. 

Table 1 . Paired t-test for difference of Leksell GammaPlan and MVEMT 



All pooled data(n=500) 



Leksell GammaPlan (mean±SD) 
26.01918.560 Gy 



MVEMT system (mean±SD) p value 

26.07518.808 Gy 0.344 



MVEMT : modified variable ellipsoid modeling technique, SD : standard deviation 



Table 2. Maximal dose and its location 



Patient 


Maximal dose (Gy)* 
(location by MVEMT) 


Maximal dose (Gy)* 
(location by LGP) 


D (max, MVEMT) - 
D (max, LGP) (%) 


As (mm) 


1 


39.94 (85.6, 108.8, 73) 


39.99 (85.6, 108.8, 73) 


-0.1 


0 


2 


33.98(101,99.6,91.6) 


32.00(101,99.6,91.6) 


6.1 


0 


3 


40.87(109.5,104.4, 63.9) 


39.93(109.5, 103.9,63.9) 


2.3 


0.5 


4 


40.95 (90, 86.6, 120.1) 


39.88 (90, 85.4, 120.1) 


2.6 


1.2 


5 


32.47(75.7, 125.4,51.9) 


31.97 (75.7, 123.8,51.9) 


1.5 


1.6 


6 


39.01 (142.5, 151.1,84) 


40.00 (142.5, 151.1, 84) 


-2.4 


0 


7 


27.00 (84.7,91.7, 110.5) 


27.00(92.7,92.7, 119.5) 


0.0 


12.08 


8 


46.34(99.7,71.9, 55.6) 


46.00 (99.7, 71.9, 55.6) 


0.7 


0 


9 


39.30 (120, 82.8, 88.5) 


40.00(120, 82.8, 88.5) 


-1.7 


0 


10 


39.82(100.1,81,58.1) 


40.00(100.6, 80.5, 58.1) 


-0.4 


0.7 


lean±SD 


37.97±5.42 


37.68±5.57 


1.78±1.80 


1.61±3.72 



*p=0.323. MVEMT : modified variable ellipsoid modeling technique, LGP : Leksell GammaPlan, SD : standard deviation 



deviation of 1.0987, a lower limit (2SD) was -2.293 with a upper 
limit was 2.014 (Fig. 2). Table 1 shows results of the paired t-test 
performed to evaluate for the differences between the LGP and 
MVEMT for all 500 points. No statistical difference was ob- 
served. 

The maximal doses according to the MVEMT were averaged 
as 37.97±5.42 Gy, while that determined by the LGP was 37.68± 
5.57 Gy (Table 2). Differences in maximal doses ranged from 
-2.4% to 6.1% (Table 2). Close agreement was observed when 
comparing both methods with an average difference of maxi- 
mal dose of 1.78±1.80%. Maximal dose points were identical 
for five cases including case 6 (Table 2). Data from case 7 showed 



largest distance (12.08 mm) between maximal dose locations 
determined by two methods. However, there was a 2.49% dif- 
ference in the dose at the coordinate (84.7, 91.7, 110.5). The av- 
erage of distance difference between maximal dose locations 
determined by two methods for entire cases was 1.61±3.72 mm 
(Table 2). 

Target volumes ranged from 0.59 to 21.7 cc (5.41±7.40) ac- 
cording to MVEMT (Table 3). LGP showed that these volumes 
varied from 0.59 to 21.1 cc (5.38±7.27). The average volume cal- 
culated by the two methods differed by 5.1 ±5.51%. No statisti- 
cally significant difference in target volume or maximal dose 
was observed. 
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Table 3. Comparison of target volume (50%) 



Patient 


Target volume of 
MVEMT (cc)* 


Target volume of 
LGP (cc)* 


(V_mvemt-V_lgp) 
/V_lgp (%) 


1 


21.70 


21.10 


2.8 


2 


1.15 


1.30 


-11.5 


3 


0.59 


0.59 


0 


4 


8.59 


8.70 


-1.2 


5 


14.97 


15.10 


-0.8 


6 


0.60 


0.59 


1.6 


7 


3.29 


3.30 


-0.3 


8 


1.73 


1.70 


1.7 


9 


0.61 


0.60 


1.6 


10 


0.89 


0.80 


11.2 


lean±SD 


5.41±7.4 


5.38±7.27 


5.115.51 



*p=0.626. MVEMT : modified variable ellipsoid modeling technique, LGP : Leksell GammaPlan, SD : standard 
deviation 



DISCUSSION 

The actual LGP algorithm is unknown. This may pose a chal- 
lenge for designing an appropriate dose verification method for 
Gamma Knife® treatment. When a method is created, it would 
be far advantageous to obtain an algorithm with known for- 
malism with secured agreement rate with LGP. In this regard, 
previous reports on verification methods are uniformly defi- 
cient due to a lack of statistical analysis 9,11,13,14 '. Additionally, re- 
sults of these studies were compared only in under certain con- 
ditions for limited parameters although many different LGP 
parameters are instantly available in clinical settings. 

MVEMT is a method that simulates a patients head as an el- 
lipsoid. The biggest difference of MVEMT from other methods 
is parameters required for calculation are generated indepen- 
dently from LGP. For example, MVEMT measure the skull ge- 
ometry from MRI, not from the bubble head measurement. 
Also, most previously developed methods require some param- 
eters such as maximal dose point and skull measurements cal- 
culated by LGP. 

Since independent verification methods adopt different cal- 
culation methods, implicit assumptions and approximations 
that produce errors under certain conditions are inevitable 13 '. 
Thus, these modalities are not applicable to all matrix points. 
Tsai et al. 13) applied average target depth from MRI images or 
ruler measurements. Since this group used a three-dimensional 
dose profile referred in unplugged condition, plugged treatments 
may produce errors depending on the numbers and location of 
plugs 13 '. They explained that the discrepancies are due to the 
high dose gradient, especially when a plug pattern is used for the 
helmet Their method revealed a 23.1% maximal dose rate mis- 
match for certain locations. Because this method produces errors 
in off-center points, it cannot be used for multiple targets. Marcu 
et al. 11 ' used a constant radius R based on a spherical head. The 
accuracy of their method depended inversely on the deviation of 
the shape of a patients head from a sphere. The method utilized 



by Zhang et al. 14 ' adopted the tissue- 
maximum ratio for each of the 201 
beams. Their method is related to accu- 
rately to the multiple- shot treatment 
plan because the accumulated error af- 
fects the normalization factor. Overall, 
the previous methods had limitations to 
be applicable for multiple targets. 

However, a large portion of Gamma 
Knife® treatments involve multiple tar- 
gets. In the authors' clinic, around 60% 
of all brain tumor cases consist of mul- 
tiple targets, most of them were meta- 
static tumors. Previous report offered 
limited examples of verification of mul- 
tiple shot treatments 11 '. Only a few au- 
thors have reviewed results of multiple 
shot treatments. The method used by Marcu et al. 11 ' had error 
ranges of 5% for a single shot and 3% for multiple shots, and 
did not provide a description of multiple targets. Zhang et al. 14 ' 
neglected to mention whether or not their study was performed 
on multiple targets 14 '. Multiple targets are subject to greater er- 
ror because matrix points are more dispersed in the cranium. 
Minimum requirements for multiple target treatment are as fol- 
lows : 1) a dose verification method with accuracy of off center 
points; 2) the ability to calculate the dose at all points in N matri- 
ces (Nx31x31x31 points); 3) an off axis dose profiles complying 
with four collimator types; and 4) inherent error not multiplied 
over the number of targets. Unless any of these requirements 
are satisfied, dose verification cannot be performed accurately. 

In our study, MVEMT showed excellent agreement with LGP 
according to the Bland- Airman analysis. This type of analysis is 
used to compare a measurement technique related to a refer- 
ence value. The limit of agreement was considered to be two 
standard deviations within a range of 95% of the comparison 
points 3 '. This method was more appropriate than correlation 
analysis when the reference method is low in error results 7 '. LGP 
also exhibits inherent error when Gamma Knife® is used in an 
actual clinical setting. For example, bubble head measurements 
are subject to uncertainty due to orientation of the ruler against 
the surface of patients' head as well as the impact of existing 
scalp hair and its elasticity 13 '. The aim of agreement analysis is 
to verify that two methods concur such that they can be used 
interchangeably. 

According to the results of our study, MVEMT agrees suffi- 
ciently well with LGP. We selected random points from random 
matrices. The selected matrices were not chosen under certain 
constraints such as distance or shot numbers. In our investiga- 
tion, 500 points were selected randomly from 297910 points 
(10x31x31x31) dispersed in the brain of 10 patients. These 
points statistically represented all arbitrary points in the skull 
regardless of shot numbers. 
Some parameters showed greater variability For example, tar- 
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get volumes had larger differences when comparing case 2 
(11.5%), and case 10 (11.2%). The rate was exaggerated because 
it was measured in rate rather than the absolute difference in 
volume. The actual volume difference was 0.15 cc for case 2 and 
0.09 cc for case 10. Even small differences in surface contour for 
each volume can create greater differences in rate. Additionally, 
Bland- Airman plots showed better results when case 6 was not 
included. It is not clear why this particular case was associated 
with worse agreement rate. However, both calculations had the 
exact same point (distance was 0) for the maximal dose point in 
case 6. Furthermore, the difference in maximal dose rate was 
2.4% (Table 2) while the difference in target volume was 1.6% 
(Table 3). Since LGP and MVEMT generate numerous parame- 
ters, different range of agreement will be produced for each pa- 
rameter. Some parameters will have greater agreement rate 
than others. This is why we compared more parameters for 
both systems than previous studies. Additionally, this also ex- 
plains why statistical analysis is necessary when comparing two 
calculation methods. 

When we consider the limitations of verification methods, 
MVEMT provided impressive results. The treatment design ap- 
plied to all matrix points in this study was performed on multi- 
ple targets. Because MVEMT is based on single beamlet distri- 
bution, accurate calculation of dose rate can be performed 
regardless of collimator plugging. The results of this study indi- 
cate that accuracy of data from MVEMT is statistically sufficient 
for guiding any treatment design. To the best of our knowledge, 
this is the first report on clinical trial of dose verification for 
multiple target treatment with a statistical analysis. Further uti- 
lization of MVEMT is expected with comparison to newly de- 
signed verification methods. 

CONCLUSION 

In this study, the accuracy of MVEMT for evaluating multiple 
target treatments was compared to that of LGK. This method 
accurately verified important parameters of multiple target 
treatment design given by LGP, which was validated by statisti- 
cal analysis. This ensures the quality of LGK treatment for mul- 
tiple target cases. Further utilization of MVEMT as a reference 
verification method is expected for designing and/or compar- 



ing new verification methods. 
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